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Abstract. We study the system of two magnetic impurities described by a two- 
impurity Kondo model where only the first impurity couples directly to the conduction 
band, while the second impurity interacts with the first through Heisenberg exchange 
coupling in a "side-coupled" configuration. We consider various choices of the impurity 
spins (S± < S2, Si =5*2, and S± > S2) and we contrast the regimes where the inter- 
impurity exchange coupling J is either lower or higher than the Kondo temperature 
T^? of the first impurity in the absence of the second. This model is a high-spin 
generalization of the two-impurity model for side-coupled double quantum dots which 
corresponds to the simplest Si = S2 = 1/2 case, where the moments are Kondo 
screened successively in two stages for J < (the "two-stage Kondo effect"). We 
show that the two-stage Kondo screening occurs generically for S2 > Si. For Si > 1, 
the second Kondo temperature is not exponentially reduced, as for Si — 1/2, but 
is approximately a power-law function of the coupling J. Furthermore, for Si > 1 
all three scales (T^\ J, Tjp) explicitly appear in the temperature-dependence of the 
thermodynamic properties. For Si > S2, there is no second stage of screening for 
J < , however in the opposite limit J > the Kondo screening of the effective 
spin 5*i — S*2 is found. 



PACS numbers: 72.10.Fk, 72.15.Qm 
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1. Introduction 

When magnetic impurities, such as substitutional defects in the bulk or adatoms on 
the surface, couple with the conduction band electrons through an antiferromagnetic 
exchange interaction, their spin is screened in the Kondo effect and the level degeneracy 
is effectively lifted [H El El HI El El [7j . When the separation between two such impurities 
is small, the impurities interact through the RKKY interaction [8], which may lead to 
critical behaviour in some parameter regimes P, \TU\ ■ A simplified description of such 
systems is the two-impurity Kondo model jTH [121 EES EEU 13 [16]: the two magnetic 
atoms are represented by quantum spin operators which are coupled by some exchange 
interaction J, and each furthermore interacts with the conduction band electrons though 
an effective Kondo exchange coupling. With few exceptions [T71 [18] , most studies of such 
models focus on spin-1/2 impurities, while real impurities may actually have higher spin 
[3]. The same may also be the case in artificial atoms, i.e. quantum dots [19], and in 
molecules with embedded magnetic ion [20]. Due to competing interactions and the 
vastness of the parameter space, a great variety of different types of magnetic behavior 
are expected. In this work we discuss a sub-class of high-spin two-impurity models in 
which only one of the spins (Si) couples to the conduction-band, while the second spin 
(S2) is "side-coupled" to the first one. Only the Si = S 2 = 1/2 limit of this family 
has been studied so far [18], and some results are known for the case of Si = 1/2 and 
arbitrary S 2 in the related Anderson- Kondo model [21]. It is shown that the two-stage 
Kondo screening [T8l E21 [23] found in the Si = S 2 = 1/2 model is a generic feature of 
all > Si models, although for Si > 1 some qualitative differences arise. 

2. Model and method 

We consider the two-impurity Kondo model 

H = X) e * c L c fcr + Jks ■ Si + JSi ■ S 2 (1) 

ha 

Operators c\ a create conduction band electrons with momentum k, spin o G {til}) and 
energy e^, while s = {s x ,s y ,s z } is the spin density of the conduction-band electrons 
at the position of the first impurity. Operators Si = {Si >x , Si >y , Si, z } and S 2 = 
{S 2 ,x, S 2j y, S 2jZ } are the quantum-mechanical impurity spin operators. Furthermore, Jk 
is the effective Kondo exchange coupling constant, and J is the inter- impurity Heisenberg 
coupling constant. The models considered are the simplest generalization of the spin-1/2 
two- impurity models for the side-coupled impurity configuration exhibiting the two-stage 
Kondo effect [HI [221 [231 [21], which can also be found in other multi-impurity problems 
[251 26] . Two-stage Kondo screening occurs when the exchange coupling between the 
impurities is weaker than the energy scale of the Kondo screening of the directly coupled 
impurity. In such situation, the moment on the first impurity is Kondo screened at the 

Kondo temperature which corresponds to a decoupled impurity, T$ , while the moment 

(2) 

on the second is Kondo screened at some exponentially reduced temperature T^ ; . A 
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simple interpretation is that the second Kondo effect occurs due to exchange coupling 
between the side-coupled impurity and the Fermi liquid of heavy quasiparticles resulting 
from the first stage of the Kondo screening [22] . 

We solve the Hamiltonian using the numerical renormalization group (NRG) 
method [27J ESI EH]. In this approach, the continuum of the conduction-band electron 
states is discretized logarithmically with increasingly narrow intervals in the vicinity of 
the Fermi level, the problem is transformed (tridiagonalized) to the form of a tight- 
binding Hamiltonian with exponentially decreasing hopping constant which is then 
diagonalized iteratively. A very effective technique to examine the magnetic behavior of 
an impurity model consists in studying the thermodynamic properties of the model as a 
function of the temperature [27J [28l [30] ; the impurity contribution to the entropy, S imp , 
then provides information on the degeneracy of the effective spin multiplets, and the 
impurity contribution to the magnetic susceptibility, Xim P , defines the effective magnetic 
moment. 

The results reported in this work have been calculated with the discretization 
parameter A = 2 using improved discretization schemes [31], [32], without the z- 
averaging, and with the NRG truncation cutoff set at 7u>n, where un is the characteristic 
energy scale at the iV-th step of the iteration. 



3. Properties of two antiferromagnetically coupled isotropic spins 

We first briefly review some properties of a decoupled pair of magnetic impurities 
described as pure spins with spin quantum numbers Si and S 2 ; without loss of generality, 
in this section we use the convention that S 2 > Si. The results will be relevant in the 
discussion of the J — > 00 limit. 

The impurities couple antiferromagnetically into an effective S = S 2 — Si spin 
object. The inner-product of spin operators is given as 

(Si-S 2 ) = -Si(S 2 + l). (2) 

The ground state multiplet can be expressed using the Clebsch-Gordan coefficients as 

\S,M)= J2 (Si,mi,S2,m 2 \S,M) \S U m x ) ® \S 2 , m 2 ). (3) 

m\,m,2 

The projections of the spin-S* object on the constituent spin operators can be easily 
computed as the expectation values of the S Z i operators in the maximum weight states 
\S,S): 

Pi = (S,S\S Ztl \S,S) (4) 
= \ (Si,mi, S 2 , S 2 - Si - mi\S 2 - Si, S 2 - 5' 1 )| 2 mi (5) 



mi 

Si(S 2 -Si) 



(6) 



and 



Si-S 2 -V 

p 2 = (S,S\S Z}2 \S,S) (7) 
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1(5*1, rni, S 2 , S 2 - 5*1 - m 1 \S 2 - Si, S 2 - S x )\ 2 (S 2 - St- mi) 
{S 1 -S 2 ){l + S 2 ) 



(9) 



S 1 -S 2 -l 

If the Hamiltonian describing the coupling of the impurities with the host conduction 
band is of the form 

H c = JiSi • s + J 2 S 2 ■ s, (10) 

then the coupling of the effective spin takes the following form 

H^ s = J cfi S-s, (11) 

with 

Jcs = riJt + r 2 J 2 , (12) 

where = Pi/ (S 2 —Si). The ratios r, are thus the multiplicative factors which determine 
the effective Kondo exchange coupling of the composite object; they are tabulated in 
Table [TJ Note that the signs of r\ is always negative, while the sign of r 2 is always 
positive. This implies that in the side-coupled configuration discussed in this work, the 
Kondo screening of the effective spin in the J — ► oo limit occurs only if the impurity 
which couples to the conduction band is the one with larger spin; in this case the 
impurity ground-state multiplet will have spin \S 2 — Si\ — 1/2. In the opposite case, 
the exchange coupling to the conduction band is ferromagnetic and the ground-state 
multiple will have spin | S 2 — S\\. 



Table 1. Multiplicative factors which determine the effective Kondo exchange 
coupling of the composite object made of the two spins locked into a S = S2 — Si 
antiferromagnetically aligned state. 

Si S 2 n r 2 



1/2 


1 


-1/3 


4/3 


1/2 


3/2 


-1/4 


5/4 


1/2 


2 


-1/5 


6/5 


f 


3/2 


-2/3 


5/3 


f 


2 


-1/2 


3/2 


3/2 


2 


-1 


2 



4. Results 

We fix pJx = 0.1 throughout this work. The Kondo temperature of a decoupled 
impurity is thus = 1.16 x 10~ 5 iy (Wilson's definition) and it is the same for any 
value of spin Si (33J, [Ml [35] . 
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The thermodynamic properties of the system for the case when the first spin is 

51 = 1/2 are shown in Figure [fl For 5*2 = 1/2 we recover exactly the prototypical 
two-stage Kondo screening where the second Kondo temperature is given by [HJ [22] 

TP = Cl TiV C2 ^, (13) 

where c\ and C2 are some numeric constants; the scaling of with / J is shown in 
the subfigure in the bottom panel. 

For S2 > 1) we observe very similar behaviour: for J < T%\ after the initial 
screening of the first impurity, the second impurity undergoes spin-5*2 Kondo screening 
which reduces its spin by one half-unit, as if it were coupled directly to the conduction 
band. The only effect of the first impurity is thus to induce much lower effective 
bandwidth D eS oc and increased density of states 1 / p c ff ^ . This picture is 
confirmed by the scaling of with / J which is very similar for all three values of 

52 in Figure [U 



S,=l/2 S_=l S=3/2 




Figure 1. Thermodynamic properties of the isotropic two- impurity clusters with 
S\ = 1/2. We plot the impurity contribution to the magnetic susceptibility, Ximp, 
and the impurity contribution to the entropy, S lmp . The symbols indicate fits to the 
universal spin-S^ Kondo magnetic susceptibility curves, which were used to extract 
the J-dependence of the second Kondo temperature, displayed in the bottom line of 
panels. 

We now consider the case when the first spin is Si = 1. We remind the reader 
that when a single impurity with spin S > 1 couples to the conduction band, one 
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half-unit of the spin is screened in a spin-S* Kondo effect, giving rise to a spin 
5—1/2 composite object which couples with the conduction band electrons with a 
ferromagnetic effective exchange coupling, thus it remains unscreened in the ground 
state [361 EH EH EH HO, EH SI]- The situation becomes more involved in the presence 
of an additional impurity, see Figure [2] 




10~ 2 10"' 10~ 2 10"' 10 1 



Figure 2. Thermodynamic properties of the two-impurity clusters with Si = 1. For 
52 = 1 and 52 = 3/2, the scaling of the second Kondo temperature with J is shown 
in the bottom panels; it should be noted that the horizontal axis here corresponds to 
J/Tpp and that the scale is logarithmic, while in figure [1] the axis corresponded to 
the inverse T^/J, while the scale was linear. The symbols for 52 = 1/2 and 52 = 1 
corresponds to fits using the universal spin- 1/2 Kondo magnetic susceptibility, while 
the symbols for 52 = 3/2 correspond to a fit using the universal spin-1 Kondo magnetic 
susceptibility. 

For S\ = 1 and 5*2 = 1/2, in the limit J ^> T% the two spins couple into a spin- 
1/2 object which couples to the conduction band with an antiferromagnetic effective 
Kondo exchange coupling J e g = r\ J with r\ = 4/3 (see Table CD) and undergoes the 
usual spin-1/2 Kondo effect, which results in fully compensated impurity spins and a 
non-degenerate singlet ground state. The expression for J e g holds strictly only when J 
is much larger than any other scale in the problem (in particular Jk and the band-width 
W) ; for J w 60 Jk we indeed find that the Kondo temperature is Tjc = 1.7 x 10~ A W which 
agrees with the expected scale of T K w Wy/pJ(4/3) exp{-l/[pJ(4/3)]} « 2 x lfrW. 
For lower J, the effective band-width is of order J rather than W, thus the Kondo 
temperature is reduced accordingly. An example of such behaviour is shown in Figure [2] 
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for J/T£> = 100. The initial free -spin In 6 entropy entropy is reduced to In 2 at T ~ J 
upon formation of the effective composite spin. This is followed by the conventional 
spin-1/2 Kondo screening (see the fit of the magnetic susceptibility with the universal 
Kondo curves). 

In the opposite limit J <C we observe the initial spin-1 Kondo screening of the 
first spin: the magnetic susceptibility goes toward 1/2 and the entropy toward 2 In 2. 
The screening is, however, abruptly interrupted at T ~ J. This can be interpreted 
as the formation of a spin-singlet object composed from the residual spin-1/2 resulting 
from the Kondo screening and the side-coupled spin-1/2 impurity. The end result is the 
same in both large-J and small- J limits; the cross-over between the two is smooth as a 
function of J. 

For Si = 1 and S 2 = 1, the behaviour for J ^> T% is particularly simple, since 
the two spins bind at the temperature T ~ J into a singlet and they no longer play 
any role. For J <C Tj^ the Kondo screening of the Si — 1 spin into a residual spin- 
1/2 is interrupted at the temperature T ~ J. The residual spin-1/2 then binds with 
the side-coupled S 2 = 1 into a new spin-1/2 composite object. Unlike the residual 
spin-1/2 resulting from the incomplete screening of a spin-1 Kondo impurity, which 
remains uncompensated since it couples to the conduction band ferromagnetically, 
the spin-1/2 composite object that emerges in this case couples with the conduction 
band antiferromagnetically, thus at some lower temperature which we again denote 

(2) 

Tk it is compensated in a spin-1/2 Kondo effect. This thus constitutes a non-trivial 
generalization of the two-stage Kondo screening phenomenology encountered in the 

51 = 1/2 cases. The differences, however, are notable: 1) there are not two, but 
three energy scales: Tjp, where the spin-1 Kondo screening takes place, J, where this 

(2) 

screening is abruptly interrupted, and T# , where the second Kondo screening occurs; 
2) the scaling of the second Kondo temperature Tg is not exponential with 1/ J. In the 
conventional two-stage Kondo effect with Si = 1/2, the only role of the coupling J is 
to set the lower Kondo temperature; no feature is observed there in the thermodynamic 
properties of the system at T ~ J. Here, the coupling J is essential to produce a 
composite spin object which then couples antiferromagnetically with the rest of the 
system, thus this scale is directly observable as a sharp change in the effective impurity 
degrees of freedom at T ~ J. The second Kondo temperature is defined by a power-law 
with the exponent near 3, with some corrections (see the lower panels in Figure [2]). 

For 5*i = 1 and S 2 = 3/2 the results for J ^> Tjp are trivial: at the temperature 
T ~ J, the spins lock into a spin-1/2 object which couples ferromagnetically with the 
conduction band, thus the composite spin remains unscreened. This is in accord with the 
expected behavior in this limit (see section[3]). For J <C the results are, however, very 
intriguing: the Kondo screening of the Si = 1 spin is interrupted at the temperature 
T ~ J. At this point, the residual spin-1/2 couples antiferromagnetically with the 

52 = 3/2 spin into a spin-1 composite object. This composite object, interestingly, 
couples antiferromagnetically with the conduction-band electrons, which leads to Kondo 
screening of one half-unit of spin at some lower temperature which we denote, yet again, 
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as Tjp . The final residual spin- 1/2 is not compensated, since it couples ferromagnetically 
with the conduction band. Thus we again observe a two-stage Kondo effect of the same 
universality class as in the Si = S 2 = 1 case. This result may, in fact, be generalized: 
for any S 2 > Si, the impurity spins will be compensated in two screening stages (the 
compensation is only partial for S 2 ^ Si). 

To further substantiate the claim that the results are generic, we show in Figure [3] 
the results for the Si = 3/2 case. For S 2 = 1/2 and S 2 = 1, i.e. for S 2 < Si we again 
find the emergence of the Kondo screening of the rigidly antiferromagnetically bound 
S 2 — Si spin in the large- J limit, and the formation of a S = (S 2 — 1/2) — Si bound state 
at T ~ J in the small- J limit. Furthermore, in the case of S 2 > Si, the two-stage Kondo 
screening is again observed in the small- J limit, again with power-law dependence of 
T { k ] on J. 



S 2 =l/2 




icr 12 icr 9 icr 6 icr 3 

TAV 



io 12 icr 9 io~ 6 io 3 

TAV 



io 12 io- 8 

TAV 



io 12 io 9 io 6 io 3 10° 

TAV 



Figure 3. Thermodynamic properties of the two-impurity clusters with S\ = 3/2. 
For 5*2 = 1/2, the Kondo screening in the large- J limit is of the spin-1 type, while for 
52 = 1 the Kondo screening in the same limit is of the spin-1/2 type. For 52 = 3/2, 
the second stage Kondo screening in the small- J limit is of the spin-1/2 type, while 
for 5*2 = 2 it is of the spin-1 type. 



5. Discussion and conclusion 

We have shown that when a second impurity is side-coupled to a Kondo impurity with 
sufficiently small Heisenberg coupling J, the spin will be screened in two stages for all 
systems, where the spin of the side-coupled impurity S 2 is equal or greater than the 
spin of the directly coupled one Si. When Si = 1/2, the second Kondo temperature is 
exponentially reduced, while for Si > 1, it is a power-law function of the coupling J. 
The difference stems from the fact that for Si = 1/2, the second stage of the Kondo 
screening occurs with a local spin S 2 which interacts with a Fermi liquid of heavy 
electrons resulting from the first screening stage, while for Si > 1 the first screening 
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stage leaves behind a residual uncompensated spin Si — 1/2, which is an extended object. 
This residual spin then rigidly binds with the spin of the side-coupled impurity at the 
temperature scale of T « J to produce a new extended spin object which then undergoes 
Kondo screening. Similar behaviour is found in the anisotropic single-impurity Kondo 
model, where an easy-plane anisotropy leads to a formation of an extended effective 
spin-1/2 degree of freedom which is Kondo screened [35]; in this problem, the second 
Kondo temperature is a power-law function of the longitudinal magnetic anisotropy 
constant D. No theory has been devised yet to map this class of problems with 
effective extended spin degrees of freedom onto the conventional Kondo model with 
a localized spin operator, thus there is presently no analytical account of these power- 
law dependences. Nevertheless, it is clear that a power-law dependence of the second 

(2) 

Kondo temperature, i.e., T K cx exp(— l/p e s J e s) oc J a , implies an inverse logarithmic 
dependence of the effective impurity parameters, i.e., p e sJ e ff oc — 1/ln J. This form 
is suggestive of the energy dependence of the renormalized ferromagnetic exchange 
coupling of the residual spin in the underscreened Kondo model, J(uS) = l/ln(a>/To), 
where T is some low energy scale [39]. This indicates that the ferromagnetic residual 
coupling might play a decisive role in determining the total effective antiferromagnetic 
exchange coupling of the composite spin object. In this scenario, the bare parameter 
J leads to the emergence of the composite spin object by antiferromagnetic binding 
of the residual spin with the side-coupled impurity spin, which occurs on the energy 
scale of uj — J, while the coupling of this object with the surrounding electron liquid is 
controlled solely by J{uj = J) = 1/ ln( J/T ). This may be explained by the fact that 
the side-coupled impurity interacts with the electron liquid only indirectly through the 
first impurity and that the sign of the relevant multiplicative factor ri is always negative 
(see Tabled]), thus the sign of the effective exchange interaction is flipped. 

In conclusion it may also be remarked that a common feature of all models 
considered in this work is that the ground state in no way depends on the J/T^ 
ratio; for any non-zero Heisenberg coupling between the impurities we always end up in 
the same fixed point, only the temperature dependence of the spin-compensation differs 
greatly. This no longer holds for problems with additional magnetic anisotropy terms 
(i.e., two-impurity extensions of models studied in Refs. [35] H2]), where level crossings 
may also occur as a function of J. This behavior will be addressed in future works. 
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